Planar resonant periodic orbits in Kuiper belt 

dynamics 

George Voyatzis and Thomas Kotoulas 
University of Thessaloniki, Department of Physics 
GR-541 24 Thessaloniki, Greece 
email: voyatzis@auth.gr, tkoto@skiathos.physics.auth.gr 

February 2, 2008 

Abstract 

In the framework of the planar restricted three body problem we study a con- 
siderable number of resonances associated to the Kuiper Belt dynamics and located 
between 30 and 48 a.u. Our study is based on the computation of resonant periodic 
orbits and their stability. Stable periodic orbits are surrounded by regular libra- 
tions in phase space and in such domains the capture of trans-Neptunian object is 
possible. All the periodic orbits found are symmetric and there is evidence for the 
existence of asymmetric ones only in few cases. In the present work first, second 
and third order resonances are under consideration. In the planar circular case we 
found that most of the periodic orbits are stable. The families of periodic orbits are 
temporarily interrupted by collisions but they continue up to relatively large values 
of the Jacobi constant and highly eccentric regular motion exists for all cases. In 
the elliptic problem and for a particular eccentricity value of the primary bodies 
the periodic orbits are isolated. The corresponding families, where they belong to, 
bifurcate from specific periodic orbits of the circular problem and seem to continue 
up to the rectilinear problem. Both stable and unstable orbits are obtained for each 
case. In the elliptic problem the unstable orbits found are associated with narrow 
chaotic domains in phase space. The evolution of the orbits, which are located in 
such chaotic domains, seems to be practically regular and bounded for long time 
intervals. 
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1 Introduction 



The dynamics of the long term evolution of small bodies revolving beyond the Neptune's 
orbit and forming the Edgeworth-Kuiper belt is a hot research topic in planetary science. 
Since 1992, when the first small body was discovered (Jewitt and Luu, 1993), hundreds of 
such small bodies, called trans-Neptunian objects (TNOs), have been detected and are in- 
cluded in the associated list of the Minor Planet Center ( |http:/ /cfa-www.harvard.edu/iau/ 
mpc.html, Jewitt, 1999). The study of their dynamical evolution provides valuable infor- 
mation towards the understanding of the origin, the formation and the structure of our 
Solar system (Morbidelli et al, 2003). The first theoretical and numerical studies (e.g. 
Torbett and Smoluchovski, 1990; Knezevic et al, 1991; Duncan et al., 1995; Gallardo and 
Ferraz-Mello, 1998; see also the review paper by Morbidelli, 1999, and references therein) 
showed the important role of resonances for the capture and the long term stability of 
TNO's. 

The most simple but efficient model for understanding the underlying dynamics of 
the resonant orbits in Kuiper Belt is the restricted three-body problem (RTBP) where 
the primary bodies are the Sun and Neptune and the small body (of negligible mass) 
represents a potential TNO. The dynamics can be understood by studying the phase 
space structure of this simple model, where chaos and order coexist. The simplest model, 
namely the planar circular RTBP, is a conservative system of two degrees of freedom and 
for each energy level the particular qualitative characteristics of the dynamics are revealed 
clearly by calculating the corresponding Poincare surface of section. Based on this tool, 
Malhotra (1996) studied some main mean motion resonances with Neptune and showed 
that resonances are associated with stable motion which is a libration with respect to the 
corresponding resonant angle a = pX' — qX — {p — q)w, where p/q denotes the external 
resonance {p < q). This property had been also indicated, using a different approach to 
the problem, by Morbidelli et al, 1995. The center of libration is a fixed point on the 
Poincare section that corresponds to a resonant stable periodic orbit for the particular 
energy level. It is well known that the periodic orbits and their stability are essentially 
associated with the structure of the phase space (Berry, 1978; Hadjidemetriou, 1993, 1998). 
Stable periodic orbits are surrounded by invariant tori and the motion is regular while 
unstable periodic orbits are associated with a hyperbolic structure. Consequently the 
usefulness of the computation of families of periodic orbits is to understand the dynamics 
in an open interval of energy levels without spending an exaggerative computational time. 
In the planar elliptic problem, which is of two degrees of freedom but non-autonomous, 
the Poincare sections are four dimensional and, subsequently, their usefulness is reduced 
significantly compared with that of the autonomous case. But the knowledge of periodic 
orbits provides essential information for the dynamics as in the autonomous case. It is 
clarified in Section 2 that such periodic orbits are isolated in a hyperplane of constant 
eccentricity of Neptune and, thus, information is provided only for a small subset of 
the phase space. Nevertheless, the localization of a stable periodic orbit in this case is 
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associated with the existence of regular orbits in the neighborhood of such a periodic orbit. 
Additionally, the numerical simulations show that the eccentricity of Neptune introduce 
small perturbations to the circular problem that do not affect significantly the long term 
stability of the resonant librating motion near the stable periodic orbits of the circular 
problem (Malhotra, 1996; Kotoulas and Voyatzis, 2004a). 

The periodic orbits of the planar circular RTBP have been extensively studied since 
the decade of '60. A detailed review, which includes the associated definitions, theoretical 
aspects, references and examples, is given by Henon (1997). The basic theory and the nu- 
merical methods for the computation of periodic orbits in the planar elliptic problem have 
been described in detail by Broucke (1969a,b) and their application in our Solar system 
and the asteroidal motion has been shown (e.g. Hadjidemetriou, 1988,1992). Applications 
also can be found for the dynamics of extra-solar systems (Haghighipour et al, 2003). As 
far as the Kuiper belt dynamics is concerned, the families of symmetric periodic orbits in 
the 1/2, 2/3 and 3/4 resonances have been studied (Kotoulas and Hadjidemetriou, 2002) 
and, as well, the 1/2, 1/3 and 1/4 resonant symmetric and asymmetric families in the 
planar circular problem (Voyatzis et al 2004). Also, the 2/3 resonant periodic orbits are 
studied by Varadi (1999), for various mass ratio and eccentricity values of the primaries. 

The aim of this paper is to study the resonant motion through the computation of 
the periodic orbits. A considerable number of resonances that are present in the Kuiper 
belt dynamics and located between 30 and 48 a.u. is included in our study. Namely, we 
examine the first order resonances p/q = m/{m + 1), where m = 1,2,..., 7, the second 
order resonances 3/5, 5/7 and 7/9 and the third order ones 4/7, 5/8 and 7/10. At this 
stage we restrict our study to the planar circular and planar eUiptic case of the RTBP. In 
section 2 we present the formulation of the RTBP model and the periodicity conditions 
used in the present work. In section 3 we present the families of periodic orbits computed 
for the planar circular model. In section 4 we discuss on the continuation of periodic orbits 
in the elliptic case and present the results of our computations. Finally, we summarize 
our results and discuss about the resonant structures obtained and their consequences in 
the long term stabihty of the small bodies. 

2 System configuration and periodicity conditions 

We consider the planar restricted three body problem with primaries the Sun and Neptune 
of mass 1 — IJ> and /j, respectively. The gravitational constant is set equal to unit. In the 
inertial orthogonal frame OXY, the motion of Neptune is either circular or elliptic round 
the center O with eccentricity e', semimajor axis a' = 1 and period T' = 27i. In the 
rotating orthogonal frame of reference Oxy, where the Sun and Neptune define the Ox- 
axis, the motion of the small body is described by the Lagrangian (Roy, 1982) 

L = lix' + y' + {x' + y')9' + 2{xy - xy)9) + ^ + ^, (1) 
2 ri ro 
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where rf = (x + /xr)^ + y'^, = (x — 1 + /xr)^ + y^, r is the distance between the primaries 



For the mass of Neptune we use the value /i=5.178xl0~^. 

In the circular problem a symmetric periodic orbit can be defined by the initial con- 
ditions x(0) = xo, y{0) = 0, i;(0) = 0, and y{0) = yo and the periodicity conditions 
are 



where T is the period of the orbit which is unknown. We can represent such a periodic orbit 
as a point in the plane Xq — h, where h is the corresponding Jacobi constant. By varying 
the value of xq (or h) we get a monoparametric family of periodic orbits with parameter 
Xq (or h). Generally, the period T changes along the family and the eccentricity cq that 
corresponds to the initial conditions of the periodic orbit is either almost constant and 
approximately equal to zero or it changes significantly along the family (see Section 2; 
Henon, 1997). The multiplicity of a periodic orbit is defined as the number of crosses of 
the orbit with the axis y = in the same direction (y > or y < 0) in a period. In the 
elliptic case, the periodic orbits are isolated points in the plane xq — yo for a particular 
value of e' and their period is T = 2/c7r, k = 1,2,... Therefore the same periodicity 
conditions (jH)) hold but the value of T is known a priori. A monoparametric family of 
symmetric periodic orbits is formed by analytic continuation varying e' (the parameter of 
the family). Such a family can be represented as a curve in the 3D space xq — yo — e'. 

In our computations the periodicity conditions are solved by a Newton- Raphson shoot- 
ing algorithm (Press et al. 1992) with accuracy 10~^^ or 10~^^ for the circular and the 
elliptic case respectively. Only for some difficult cases (e.g. near collisions) we were forced 
to decrease the accuracy by one decimal order. The numerical integrations were performed 
using a control step Bulirch-Stoer method. The starting points for the computation of 
the families of periodic orbits in the circular case are obtained from Poincare sections 
at a particular level of the Jacobi constant. For the elliptic problem we started from the 
known bifurcation points derived from the families of the circular problem. Our algorithm 
controls automatically the possible change of the multiplicity of the orbits and the step 
of the parameter of the family. Since a symmetric periodic orbit crosses the y = axis at 
least twice, the algorithm may change the starting point of the periodic orbits when the 
convergence to the prescribed accuracy fails. The linear stability of the periodic orbits 
is determined from the corresponding stability indices. For the planar circular case the 
stability index k = an + 022 is computed, where aij are the elements of the monodromy 
matrix (Henon, 1997). For the planar elliptic case we calculate the indices ki and k2 intro- 
duced by Broucke (1969a,b). For both cases the computation is based on the numerical 
solution of the variational equations. 



and 6 the angle between the axes Ox and OX. For the circular problem it is r = 1, 6* = 1 
and there exist the Jacobi constant h written in the form 




y{0) = y(T/2) = 0, x(0) = x(T/2) = 0, 



(3) 
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Figure 1: Characteristic curves showing the continuation of famihes of elhptic periodic 
orbits of first order from the circular ones (/i 7^ 0). The almost straight line segments, 
indicated by the circles, consist the family of circular orbits (first kind). 

3 The planar circular problem 

We consider the planar circular case of the RTBP. The families of symmetric periodic 
orbits are classified in two different kinds (Poincare, 1892; Henon, 1997): 
Families of circular orbits (or first kind) : The periodic orbits correspond to nearly circular 
orbits for the small body. The period T, and subsequently the resonance n/n' = T/T', 
varies along the family. 

Resonant Families (or second kind) : The periodic orbits correspond to almost elliptic 
orbits for the small body. The eccentricity Cq increases along the family but the ratio n/n' 
is almost constant and rational n/n' ^ p/q, p,q ^ Z. 

In the exterior resonances n/n' = p/q, studied in this paper, it is g > p and the 
difference A = q — p defines the order of the resonance. The circular periodic orbits of the 
unperturbed problem (/i = 0) are continued for yU > and a family C of first kind exists 
at each particular value of fi. The families of resonant periodic orbits originate from the 
p/g-resonant orbits of the family C. In each resonance, independently of its order, there 
exist two different families, / and II of second kind, which differ in phase. On family / 
the small body is initially at perihelion and on family II it is at aphelion. The families 




5 




Figure 2: Families of resonant periodic orbits for the ) 4/5 and b) 5/6. The symbols 

X denote close encounters with Neptune while the encircled ones denote close encounters 
with the Sun. The segments ll^^ and J^yg are the only that consist of unstable periodic 
orbits and are presented by a thin curve. 

may consist of many segments separated by gaps which correspond to close encounter 
orbits where computations fail. We use the notation I^^^ (or Hp/q), where p/q denotes 
the corresponding resonance and n numbers the possible different segments in the family. 
In the following discussion and plots we select the variable Xq as the parameter for the 
characteristic curves that present the resonant families. 

3.1 1st order resonances 

We consider the first order resonances with p/q = m/{m + 1), m = 2,3, ..,7. In this 
case the circular family of the unperturbed problem breaks close to the resonances when 
fj, > and the persisting segments are separated by gaps. These segments, consisting 
of almost circular orbits, are continued smoothly giving rise to the resonant families / 
and //. The six characteristic curves computed are presented in Fig. 1. We obtain that 
a family /m/(m+i) joins smoothly (through a segment of almost circular orbits) with the 
family //(m+i)/(m+2)- By decreasing the parameter xq, the famihes / approach a collision 
orbit with Neptune for h ~ —1.5. The first order resonances are considered exceptional 
because they follow the above scenario instead of a regular bifurcation (Guillaume, 1974; 
Hadjidemetriou, 1993). In our case a rather different structure is observed for the families 
/g/r and IIj/s because the resonant family Iq/j avoids the collision with Neptune. A 
new resonant family of unstable periodic orbits, denoted by Hys, is found and the three 
families and Hyg) form a closed characteristic curve. 

In Fig. 2a and 2b we present the variation of the resonance and eccentricity, re- 
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Figure 3: The variation of the resonance (left) and the initial eccentricity of the periodic 
orbits along the families (the parameter of the family is considered to be the variable x). 



spectively, along the first order resonant families shown in Fig. 1. Actually we present 
the mean values, computed along one period time interval. We observe that the ratio 
n/n' vary along the part of the family where the orbits are almost circular. The length of 
these parts, which corresponds to the the plateau e ~ of the curves in Fig. 2b, decreases 
rapidly as m increases, and disappears for the closed characteristic curve mentioned above. 

In Fig. 1 only the first segment of the families and II^/{m+i) shown in 

order to emphasize their origin. In Fig. 3 the complete resonant families 4/5 and 5/6 are 
presented. For all the resonances with 1 < m < 6 the family segment Im/{m+i) consists 
of unstable periodic orbits of multiplicity one and terminates at a collision orbit with 
Neptune. The family continues with new segments separated by collisions with Neptune, 
and extends up to a collision with the Sun. For m < 4 (i.e in the resonances 2/3 and 3/4) 
only one collision orbit is obtained while for m > 4 two collision orbits exist. In all cases, 
the segments after the first collision consist of stable periodic orbits. In the families II, 
the segments IIm/{m+i) < m < 7) reach a collision orbit too, as h increases. Beyond 
this point the family continues with the segment IIm/(m+i) "which either extends up to a 
collision orbit with the Sun (cases for 1 < m < 4) or is interrupted by a second collision 
orbit with Neptune (cases 5/6 and 6/7). All the families IIm/{m+i) consist of stable orbits. 
Only the periodic orbits which are close encounter orbits prove to be unstable but this 
result may be an artifact of the limited accuracy of computations. 
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Figure 4: Families of periodic orbits for the second order resonances. Bold or thin curves 
denote stable or unstable orbits respectively. The family C of circular orbits is also 
indicated. The symbols x denote close encounters with Neptune while the encircled ones 
denote the termination of the family to a collision orbit with the Sun. 

3.2 Second order resonances 

We consider the second order mean motion resonances p/ q=3/5, 5/7 and 7/9. The families 
of resonant periodic orbits bifurcate from the circular family C by analytical continuation 
for n ^ 0. The bifurcation points are those on family C where n/n' = p/q. For each 
resonance two families bifurcate, called family / and family //. The periodic orbits 
included in family / cross vertically the axis Ox (i.e. xq = 0) only for x > 0, while in 
family // such a vertical cross occurs only for a; < 0. Thus, the two families are presented 
in separated plots Xq — h shown in Fig. 4. 




(a) 






Figure 5: Periodic orbits of the family I^/^ in the x — y rotating frame, a) Xq = 1.5, 
b) Xq = 1.77 (before the collision with Neptune) c) Xq = 2.0 (after the collision) and d) 
Xq = 2.7 (close to collision with the Sun). The position of the Sun (x = /i) and Neptune 
(x = 1 — /x) is indicated. 
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The resonant family Jp/g starts from the bifurcation point having unstable periodic 
orbits and is interrupted by a close encounter with Neptune for h —1.5. After a close 
encounter the families continue with stable orbits. In the 3/5 resonance, this continuation 
extends smoothly up to high eccentricity values (i.e. up to a collision orbit with the Sun). 
For the 5/7 resonant case, the family reaches to a second close encounter with Neptune 
and then terminates at a collision orbit with the Sun. In the 7/9 resonant case three 
collision orbits occur along the family. The multiplicity of the orbits is affected when a 
close encounter takes place. An example of some typical 3/5 resonant orbits of the family 
/ is shown in Fig. 5. 

Excluding close encounter orbits, family Hp/q consists of stable periodic orbits. Along 
the families we obtain one, two and three close encounters with Neptune for the 3/5, 5/7 
and 7/9 resonances respectively. All families terminate at a collision orbit with the Sun 
(Fig. 4b). 




Xq 



Figure 6: a) Families of third order resonances 4/7, 5/8 and 7/10. The symbol x and the 
encircled x indicate collisions with Neptune and the Sun respectively. The solid circles 
indicate almost double collisions, b) A stable periodic orbit of the family //7/10 close to 
a double collision {h —0.75) in the rotating x — y frame, c) The same orbit as in (b) 
presented in the inertial frame X-Y. The circular orbit N of Neptune is also plotted. 
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3.3 Third order resonances 



Our study includes the resonant cases p/g=4/7, 5/8 and 7/10. The dynamics of third 
order resonances shows similar qualitative characteristics to these of the second order. The 
families / and // of elliptic resonant periodic orbits bifurcate from the circular family and 
they are presented in Fig. 6a. Families / start with unstable orbits and extend up to 
close encounter with Neptune at ~ —1.50. The families continue after this point with 
stable orbits and extend up to a second close encounter. At this point the family 74/7 
seems to terminate or become strongly chaotic and is difficult to be localized by the 
computations. This is also the case for the family /7/10 after the third close encounter 
while the continuation of the family /s/g is possible up to a collision orbit with the Sun. 

The periodic orbits of families Hp/q are all stable (except in the neighborhood of 
collisions). Along the families //s/g and //7/10 two close encounters with Neptune occur 
but only one for the family 1/4/7. The family 1/5/8 is continued up to very high eccentricity 
values and terminates at a collision with the Sun. The families //4/7 and //7/10 extend up 
to high eccentricity values and their computation terminates when close encounters with 
both Neptune and Sun take place. A linearly stable periodic orbit of the family //7/10 
that approaches the above double collision is shown in Fig.6b and 6c in the rotating and 
inertial frame respectively. The intersection of the orbits, which is shown in Fig. 6c, is a 
usual feature of resonant motion. However, because of resonance, there is phase protection 
and coUisions are avoided (Hadjidemetriou 1988, Morbidelli, 1999). 



4 The elliptic restricted three-body problem 

The resonant periodic orbits of the circular problem are bifurcation points (BPs) or, in 
different terminology, generating orbits for the elliptic problem (e' 7^ 0) when their period 
is multiple of T' = 27r. Namely, starting from a generating orbit at e' = 0, a family of 
periodic orbits of constant period is formed by varying e' (the parameter of the family). 
In Table 1 wc present the values of eccentricity and Jacobi constant that correspond to 
the bifurcation points found in the studied resonances (see also Kotoulas and Voyatzis, 
2004b). From each BP/ two families bifurcate, denoted as E^^"^ or Ef^'', where p/q is the 
resonance, / is the bifurcation point (according to the numbering in the first row of Table 1, 
which is based on the ascending sorting of the corresponding eccentricity values) and the 
second subscript p or a denotes the initial position of Neptune is at perihelion or aphelion 
respectively. The case n = is discussed in section 4.2. In Table 1, the symbol "S" or 
"U" denotes the stability of the bifurcating periodic orbits at e' ^ (stable or unstable 
respectively). The first symbol refers to family Ef^*^ and the second one to family E^^''. 
In ail cases, the indicated stability is preserved along the family up to e' = 0.01, which 
is the actual eccentricity of Neptune. But we calculate the families and their stability up 
to high values of e' in order to provide some general results on the resonant structures in 
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the elliptic RTBP. In this case Neptune refers to a fictitious planet. 



Table 1: The bifurcation points (BP/) where the families E^^'^ and E^l'^ of periodic orbits 
originate. We present the corresponding initial eccentricity value Cq and the Jacobi con- 
stant h in brackets. The symbol "S" or "U" denotes the stability character of the periodic 
orbits close to the bifurcation point. 



n/n' 


BPO 


BPl 


BP2 


BP3 


BP4 


2/3 


_ 


0.469 SU 


_ 


_ 


_ 




_ 


(-1.393) 


_ 


_ 


_ 


3/4 


_ 


0.329 SU 


_ 


_ 


_ 




_ 


(-1.452) 


_ 


_ 


_ 


4/5 


_ 


0.253 SU 


0.871 SU 


_ 


_ 




_ 


(-1.473) 


(-0.960) 


_ 


_ 


5/6 


_ 


0.205 SU 


0.749 SU 


_ 


_ 




_ 


(-1.483) 


(-1.146) 


_ 


_ 


6/7 


_ 


0.172 UU 


0.649 SU 


0.960 SU 


_ 




- 


(-1.488) 


(-1.253) 


(-0.743) 


- 


3/5 


0.0 uu 


0.427 US 


0.800 SU 








(-1.541) 


(-1.428) 


(-1.065) 






5/7 


0.0 uu 


0.278 SU 


0.562 SU 


0.778 US 


0.936 SU 




(-1.518) 


(-1.474) 


(-1.325) 


(-1.102) 


(-0.792) 


7/9 


0.0 uu 


0.203 SU 


0.427 SU 


0.606 US 


0.766 SU 




(-1.510) 


(-1.487) 


(-1.406) 


(-1.287) 


(-1.122) 


4/7 


0.0 UU 


0.027 UU 


0.400 UU 


0.900 SU 






(-1.549) 


(-1.549) 


(-1.408) 


(-0.864) 




5/8 


0.0 UU 


0.029 SU 


0.335 US 


0.800 SU 






(-1.535) 


(-1.535) 


(-1.468) 


(-1.067) 




7/10 


0.0 UU 


0.025 SU 


0.249 US 


0.905 SU 






(-1.520) 


(-1.520) 


(-1.485) 


(-0.872) 





4.1 First order resonances 

The computations reveal that there exist bifurcation points (BPs) of periodic orbits in the 
elliptic problem, only in the families Ilp/q, {p/q = ni/ (m + 1). rn = 1, 6). Each one of 
the 2/3 and 3/4 resonances has one BP and two families originate from each one. These 
families were studied in Kotoulas and Hadjidemetriou (2002). In the 4/5 resonance there 
exist two bifurcation points. The families E^^^ and E^^J' start from BPl and the Et,^ and 
E2a start from BP2. The families E\^ and E2^ include stable orbits, while the orbits of 
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(a) 



Figure 7: a) The resonant families, with e' as a parameter, of the elhptic problem at 
the resonance 4/5. The families bifurcate from the family 1/4/5 located at e' = and 
presented by the thin curve, b) The stable periodic orbit of the family Ef^^ for e' = 0.4 in 
the rotating frame Oxy. c) the same as in (b) for the family -E^^^. This orbit is unstable. 

the families Ef[^ and -^20^ unstable. These families are presented in the Xq — i/Q — e' 
space in Fig. 7a. In the same figure the corresponding family of the circular problem is 
also indicated in the plane e' = 0. Two samples of periodic orbits of the families E\^^ and 

EZ' for e' = 0.4 are shown in Fig. 7b, c. Without loss of important information, it is more 
convenient to present the families in the projection space e' — Cq, where Cq = e(0) is the 
eccentricity value of the orbit that corresponds to the initial conditions of the periodic 
orbit. Certainly, a point in the plane e' — Cq does not represent a unique orbit. The 
variation of the eccentricity e = e(t) along a periodic orbit is rather insignificant and 
evidently periodic. Therefore, we consider the value cq as the eccentricity of the whole 
periodic orbit. By using the above mentioned plane, the computed families are shown in 
Fig. 8. Sohd and dashed curves refer to the families E^j^'^ and E^l'^ respectively. Also 
the type of stability is indicated. Thin lines denote stable orbits and bold ones denote 
unstable orbits. 

In Fig. 8 we observe that in all resonances and along the families Ena , n = 1,2, .., 
the initial eccentricity (eo) increases, as e' increases, monotonically. The families seem 
to continue up to the rectilinear case (e' = 1). But in most cases and for e' > 0.9 
our computations fail to localize the periodic orbits for various reasons (e.g. the y or 
its derivatives with respect to y or a; become large). Exceptional cases are the families 
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Figure 8: The presentation of the famihes E'^^ (sohd curves) and i?;*^ (dashed curves) of 
first order resonances in the plane e' — eo, where cq = e(0). Bold curves indicate the family 
segments of stable periodic orbits and thin curves indicates the segments of unstable ones. 



-^2a 5 -^2a i?3^ , which start from the BP with the largest eccentricity value. These 
families tend to terminate at a collision orbit with the Sun as e' increases. For the families 
'T- = 1,2, .. we observe that they start with a decreasing cq as e' increases and the 
orbits are almost elliptic with longitude of pericenter u; = 180°. If the orbits attain the 
value Co = 0, at some particular value e', the family continues with orbits of u) = 0° and 
the initial eccentricity increases along the rest of the family. In all families, the stability 
type is preserved except in families E'^^ and E^^,^ . Especially, in family E\^^ we observe 
a segment of unstable orbits in the interval 0.17 < e' < 0.55. The end points of such 
an interval are candidate bifurcation points for families of asymmetric orbits (Beauge et 
al,2003; Voyatzis et al 2004). 
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Figure 9: Families of periodic orbits of the 3/5 and 5/7 resonances. The presentation is 
the same as that in Fig. 8. 



4.2 Second order resonances 

In the second order resonances (p/g, q = p + 2) both famihes I and // of the circular prob- 
lem have bifurcation points, which are starting points for families of the elliptic problem. 
Additionally, bifurcation points (denoted as BPO in table 1) exist in the circular family 
C where the period is T = gvr. Such circular orbits are continued for e' 0, providing 
periodic orbits of period T = 2qTT (Hadjidemetriou 1992). There are two different initial 



configurations giving rise to these families denoted as Eq['' and 



Family E^i" 



02 



Family E^^'^ 



Sm-Nper-Bper (t 



0) ^ SVN-Nap-Bap {t = T/2) 

02 . Bap-SVN-Np,r (t = 0) ^ Sp,,-SUN-iV,p {t = T/2) 
where B and denote the small body and Neptune respectively and the subscripts ap 
and per denote the position of the bodies (at apocenter or pericenter respectively). 

In the case of the 3/5 mean motion resonance there exist three bifurcation points (see 



Table 1). From BPO, which belongs to the family C, the families Eq(^ and -E'02"' bifurcate 
with unstable orbits. The BPl belongs to the family 7/3/5 of the circular problem and 
BP2 belongs to a stable segment of the 13/5 one. In Fig. 9a the corresponding families 
are presented in the e' 



'3/5 

Co plane and their stability is indicated. 



As in the case of 



first order resonances, along the families Ef^ or E^^ that bifurcate from the family 
II the eccentricity eo increases or decreases, respectively, as e' increases. The opposite 
situation holds for the families that bifurcate from the / family of the circular problem. 
In comparison to first order resonances, a different structure is observed for the families 
Eq2^ and -Eip^, which join smoothly at e' ~ 0.85 and can be considered as one family 
starting from BPO and ending at BPl. 
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In the case of 5/7 resonance, there exist five bifurcation points. The bifurcation point 
BPO belongs to an unstable segment of the family C of the circular problem. The points 
BPl and BPS belong to the family 1/5/7 and the BP2 and BP4 belong to the Js/y one. 
The corresponding families and their stability is shown in Fig. 9b. As in the case of 3/5 
resonance, the families E02 and Eip join smoothly. The type of stability of orbits changes 

5 /7 

twice along the family E'gp a segment of unstable orbits exist. The edges of this 

segment may be bifurcation points for families of asymmetric orbits, as mentioned above 
for the E^p^ family. 

The 7/9 resonance shows five bifurcation points distributed as in the 5/7 resonant case. 
Also, the generated families have the same qualitative features as the ones presented for 
the 5/7 resonance. 
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Figure 10: a) Families of periodic orbits of the 4/7 resonance in the planar elliptic RTBP. 
b) The evolution of the eccentricity along some orbits which start near the periodic orbits 
Tri (unstable), Tr2 (stable) and Tr^ (doubly unstable) indicated in (a). Note the different 
scales in the vertical axis. 



4.3 Third order resonances 

The families of the circular problem at the resonances p/q =4/7, 5/8 and 7/10 exhibit four 
points from which families of periodic orbits of the elliptic problem bifurcate (see Table 
1). In all cases, the bifurcation point BPO belongs to the family C and corresponds to a 
periodic orbit of period 2g7r/3. This orbit is continued to the elliptic problem if we assume 
that it is described three times providing periodic orbits of period T = 2g7r. The rest of 
bifurcation points (BPl, BP2 and BPS) belong to the family //. The bifurcating families 
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for the 4/7 resonance and their stabihty type are shown in Fig. 10a. The famihes E^g^ 
and EqI'^ join smoothly at a relatively low value of the primary's eccentricity (e' ~ 0.13). 
Such a property is also obtained for the famihes Efp"^ and £'2^^, which join at e' f« 0.8. 

However, the family E2U does not show a typical continuation. It starts out from BP2 
with unstable orbits and continues by increasing e' up to the value 0.767. At this point, 
e' decreases along the family and the periodic orbits arc stable. The stable orbits occupy 
a short segment (0.745 < e' < 0.767) in the family which is followed by a segment with 
complex unstable orbits (0.740 < e' < 0.745). For e' < 0.74 the orbits become doubly 
unstable and the family terminates at e' ~ 0.296. This termination point is not either a 
collision orbit or a strongly chaotic one and we are not able to interpret it. Also, there is 
not any obvious reason for the failure of calculations. 

Comparatively to the 4/7 resonance, the resonances 5/8 and 7/10 have four bifurcation 
points too. However, the stability type of the bifurcating families differs and is indicated 
in Table 1. The continuation of the families shows the same quahtative characteristics as 
these described above. 

Closing our study on the elliptic problem we should state that in all cases of unstable 
periodic orbits the stability indices have values slightly larger than 2.0. This fact, beside 
numerical integrations, suggests that chaos is not present in a sense of practical impor- 
tance. Nevertheless, the presence of stable and unstable asymptotic manifolds changes 
the phase space topology. The orbits that start in the neighborhood of unstable periodic 
orbits show a large variation in their orbital elements in comparison with that shown by 
the orbits which start near a stable periodic orbit. A typical example of the evolution 
of the eccentricity along some orbits of the same family E2U but of different stability 
type is shown in Fig. 10b. The presented trajectories correspond to the initial conditions 
xo,yo — and yo which coincide with that of the periodic orbits Tri, Tr2 and Trs in- 
dicated in Fig. 10a while we have set Xq = 10~^. The time scale is based on Neptune's 
period, namely 27rt.u.=165 years. For the trajectory starting near Tr2, which is stable, 
the evolution of the eccentricity e = e{t) shows small amplitude oscillations around the 
starting value. In contrast, the evolution of e{t) along the orbits that start near the simply 
unstable periodic orbit Tri and the doubly unstable Tr^ shows remarkable deviation from 
the initial value. However, for both cases of instability, the evolution seems regular. 

It is worth to mention that the orbits along the families Ef^'^ intersect the orbit of 
Neptune in general (as in the case of Fig. 6c) but the collision of the bodies is avoided 
because of phase protection. For the famihes Ef^'^ the intersection of orbits occurs only 
for orbits belonging in some particular short segments of the families. We did not find 
any colhsion orbits along the famihes studied in the elhptic RTBP. 
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5 Conclusions 



In this paper we presented families of resonant symmetric periodic orbits obtained for the 
planar RTBP and considering ^ = 5.175 • 10^"'' (the normalized mass of Neptune). Thus, 
our study on the circular RTBP is associated to the Kuiper belt dynamics and its most 
important resonances between 30 and 48 a.u. Particularly we studied systematically the 
first order resonances 2/3, 3/4, 4/5, 5/6 and 6/7, the second order ones 3/5, 5/7 and 7/9 
and the third order ones 4/7, 5/8 and 7/10 for the circular and the elliptic RTBP. Based 
on well established methods and previous results of Henon, Broucke and Hadjidemetriou, 
we determined all the families of the above mentioned resonances and presented their 
main characteristic such as bifurcation points, stability, continuation and collisions. 

For the circular planar problem we found that all studied resonances have two families 
of periodic orbits: the family / and the family II. Family // consist always of stable 
periodic orbits while unstable periodic orbits are obtained only in a short segment of family 
/ located between the origin of the family and the first collision point which is met. Also 
very short segments of unstable orbits are obtained close to collisions but in these cases the 
computations are ambiguous. The stable orbits constitute centers for resonant librating 
motion where Trans-Neptunian objects can be captured. The families, though they are 
interrupted by collisions, extend up to large values of the Jacobi constant preserving their 
stability i.e. regular highly eccentric motion can be localized for each resonance at least 
in the level of the RTBP model. For all resonances we located the bifurcation points at 
e' = which are starting points for families of periodic orbits for the elliptic problem. At 
each one of them a pair of families and Ep) originates. Thus for e' ~ 0.01, which is the 
eccentricity of the Neptune's orbit we obtain isolated periodic orbits, one in each family. 
The families of the elliptic problem start withm either stable or unstable orbits. The 
stability type of the bifurcating families is preserved up to the value e' of the Neptune's 
orbit, but in many cases it changes for larger values of e'. The corresponding eigenvalues 
of the unstable orbits are always very close to the value +1 (i.e. at least one stability index 
is slightly greater than 2). This situation has been also observed for internal resonances 
(Hadjidemetriou, 1993) and indicates the existence of weak instability at the particular 
regions in phase space. Indeed, our numerical simulations show that in the neighborhood 
of unstable periodic orbits of the elliptic problem the motion seems regular for long term 
evolution but the orbital elements show slow and large variation. Generally, collisions are 
responsible for the generation of strong chaos but in the elliptic problem all the families 
found avoid collisions with Neptune. Only the families Ea which bifurcate from the most 
eccentric bifurcation point seem to terminate at a collision with the Sun as e' — > 1. 

As far as the planar RTBP is concerned, our systematic study on resonant periodic 
orbits verifies the dynamical structures obtained in the past for other resonant cases. 
Additionally, we should remark the following points: a) The scenario of bifurcation and 
continuation of the families of first order resonances (Guillaume, 1974; Hadjidemetriou, 
1993), shows some differentiation when the resonant families avoid collisions with the 
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planet. In our study, this is obtained for the I^p and //y/s famihes, which join together 
forming a closed characteristic curve b)The famihes of the elhptic problem seem to extend 
either up to the rectilinear case (e' = 1) or join smoothly with another family forming 
characteristic curves that start from one bifurcation point and terminate to another one. 
c) It has been suggested (Voyatzis et al, 2004) that the existence of asymmetric periodic 
orbits is associated with the existence of continuous segments of unstable orbits along a 
family of symmetric orbits. This condition is true for the resonances of the form q — 
1,2,... In our study of the circular problem no any families with such a property have been 
found. In the elliptic problem the families E^^^ and include a segment of unstable 
periodic orbits and, subsequently, we may conjecture the existence of asymmetric periodic 
orbits in these cases. 

An interesting extension of the present work is to determine the families of periodic 
orbits in the three dimensional (3D) RTBP. The bifurcation points of families of 3D 
periodic orbits are given in Kotoulas and Voyatzis (2004b). The 3D resonant structure 
should provide useful information about the capture of trans-Neptunian objects in highly 
inclined orbits. Interesting structures of resonant motion can also be provided by the 
study of asymmetric resonances located beyond 48 a.u. and are not included in the 
present study. 
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